Charge order and spin-singlet pairs formation in Ti407 



I. Leonov"'^, A. N. Yaresko^, V. N. Antonov'^, U. Schwingenschlogl'*, V. Eyerf*, and V. I. Anisimov^ 
^ Theoretical Physics III, Center for Electronic Correlations and Magnetism, 
Institute for Physics, University of Augsburg, Germany 
^ Max-Planck Institute for the Physics of Complex Systems, Dresden, Germany 
^ Institute of Metal Physics, Vernadskii Street, 03142 Kiev, Ukraine 
* Theoretical Physics II, Institute for Physics, University of Augsburg, Germany and 
^ Institute of Metal Physics, Russian Academy of Science- Ural Division, 620219 Yekaterinburg GSP-170, Russia 

(Dated: February 2, 2008) 

Charge ordering in the low-temperature trichnic structure of titanium oxide (Ti407) is inves- 
tigated using the local density approximation {LT>A)-\-U method. Although the total 3d charge 
separation is rather small, an orbital order parameter defined as the difference between t2g occu- 
pancies of Ti^+ and Ti*+ cations is large and gives direct evidence for charge ordering. Ti 4s and 
4p states make a large contribution to the static "screening" of the total 3d charge difference. This 
effective charge screening leads to complete loss of the disproportionation between the charges at 
3+ and 4-1- Ti sites. The occupied t2g states of Ti^^ cations are predominantly of d^y character and 
form a spin-singlet molecular orbital via strong direct antiferromagnetic exchange coupling between 
neighboring Ti(l) and Ti(3) sites, whereas the role of superexchange is found to be negligible. 

PACS numbers: 71.20.-b, 71.28.4-d, 71.30.-|-h 
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I. INTRODUCTION 



The mixed valent transition metal oxides, that simul- 
taneously contain metal atoms in two (or more) differ- 
ent valence states, are of strong current interest^ One 
of the classical examples of such a system is magnetite 
(Fe304), in which a first-order metal-insulator transition 
occurs at ^--^120 K«2. According to Verwey, this transition is 
caused by the spatial ordering of 2-1- and 3-|- Fe cations on 
the octahedral B-sublattice of the inverted spinel struc- 
ture AB204.'^'^ Recently, a local spin density approxima- 
tion (LSDA)+C/ study of the low-temperature phase of 
Fe304 resulted in a charge and orbitally ordered insulat- 
ing ground state with a well pronounced orbital order. 
However, the strong difference in t2g occupancies of 2+ 
and 3-1- Fe was found to be drastically reduced by effec- 
tive "static" screening.^ A similar result (see Ref. ^ has 
been obtained for another iron oxide, containing both 
2-1- and 3-|- Fe cations, iron oxoborate (Fe20B03), which 
shows a broad semiconductor-to-semiconductor transi- 
tion at '^317 K associated (as in Fe304) with a spatial 
order-disorder transformation of 2-1- and 3-1- Fe cations 
on quasi-one-dimensional Fe chainsAiSiii 

The aforementioned phenomena of sharp metal- 
insulator transitions associated with pronounced charge 
and/or orbital ordering are characteristic for the Magneli 
phases Mn02n-i (M = Ti, V)ii^ These compounds form 
a homologous series and have been studied recently 
to understand the differences in crystal structures and 
electronic properties between the end members MO2 
(n — > 00) and M2O3 (n — 2)}^ In particular, the metal- 
insulator transition of VO2 discovered some fifty years 
ago still is the subject of ongoing controversy and is an- 
other "hot topic" in solid state physics. LDA calculations 
have revealed strong influence of the structural degrees of 
freedom on the electronic properties of VO2 and neigh- 



bouring rutile-type dioxides li^ji^tiSiiLi^ In this scenario 
the characteristic dimerization and antiferroelectric dis- 
placement of the metal atoms translate into orbital or- 
dering within the t2g states and a Peierls-like singlet for- 
mation between neighbouring sites. Recently, this was 
confirmed by LDA-t-DMFT calculations, which suggested 
to regard the transition of VO2 as a correlation-assisted 
Peierls-transition>i2 

Ti407 titanium oxide is another remarkable member 
of the Magneli phases with n=4 which shows metal- 
insulator transitions associated with the spatial charge 
ordering. It is a mixed valent compound which has an 
even mixture of 3-1- and 4-1- Ti cations (Ti^+Tij+Oy), 
corresponding to an average 3d occupation of 1/2 elec- 
tron per Ti site. Electrical resistivity, specific heat, 
magnetic susceptibility, and x-ray diffraction data reveal 
two first-order transitions in the temperature range of 
130-150 Ki22i2ii2a At 150 K a metal-semiconductor tran- 
sition occurs without measureable hysteresis in resistivity 
and specific heat. It is followed by a semiconductor- 
semiconductor transition at 130-140 K, which again is 
characterized by an almost two orders of magnitude 
abrupt increase in electrical resistivity and has a hys- 
teresis of several degreesi^2i2i The magnetic susceptibil- 
ity shows a sharp enhancement when heating through 
150 K. However, it is small and temperature independent 
below this temperature and does not show any anomaly 
at 140 K. 

The crystal structure of Ti407 (see Fig. can be 
viewed as rutile-type slabs of infinite extension and four 
Ti sites thickness, separated by shear planes with a 
corundum-like atomic arrangement. Below 130 K it crys- 
tallizes in a triclinic crystal structure with two formula 
units per primitive unit cell3iS4i2^ Four crystallograph- 
ically inequivalent Ti sites are found at the centers of 
distorted oxygen octahedra. They form two types of 
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FIG. 1: (color online) The low-temperature crystal structure 
of Ti407. Chains of four Ti sites run parallel to the pseudo- 
rutile c-axis. Red and blue (light and dark on the black and 
white image) chains of four Ti atoms correspond to the (a) 
and (b) chains Ti atoms, respectively. Further gradation of 
red and blue on light and dark subsets indicates inequivalent 
Ti sites in (a) and (b) chains. 



chains, namely, (a) 1-3-3-1 and (b) 2-4-4-2, which run 
parallel to the pseudo-rutile c-axis and are separated by 
the crystallographic shear planes. Although interatomic 
distances in the (b)-chain are almost uniform (3.01 and 
3.07 A between 4-4 and 2-4 Ti sites, respectively) they 
are remarkably different for the (a)-chain (3.11 and 2.79 
A between 3-3 and 1-3 Ti sites). 

Accurate determination of the crystal structure al- 
lowed to elucidate the nature of the three phases dis- 
tinguished by the two first-order transitionsi22iSi2^iS^ In 
particular, in the metallic phase the average Ti-0 bond 
lengths for crystallographically inequivalent TiOg octa- 
hedra are very similar which results in the average valence 
state of 3.5-1- per each Ti cation. Below 130 K charge 
has been transferred from the (b) to the (a)-chains. In 
addition, Ti'^"'" cations in alternate (a)-chains are paired 
to form nonmagnetic metal- metal bonds, whereas in the 
intermediate phase pairing also persists but its long- 
range order calls for a fivefold supercelliiSi Thus, the 



130-140 K transition is associated with a transition to 
the phase with a long-range order of Ti'^^-Ti'^"'' pairs, 
whereas above 150 K 3-1- and 4+ Ti cations are disor- 
dered. The presence of the Ti^+-Ti'^+ pairs strongly dif- 
ferentiates Ti407 from Fe304 and results in two steep 
first-order transitions found in the electrical resistivity. 

Recent LDA band structure calculations of both high- 
and low-temperature phase of Ti407 results in signifi- 
cant t2g charge separation between crystallographically 
independent 3-1- and 4-1- Ti sites in the low-temperature 
phase, whereas a rather isotropic occupation of the t2g 
states has been found at room-temperaturei^ While, in 
addition, an orbital order at the Ti d} chains originat- 
ing from metal-metal dimerization was found, the LDA 
gave only metallic solution with semimetallic-like band 
overlap instead of the semiconducting gap. This prob- 
lem is overcome in our work taking into account strong 
electronic correlations in Ti 3d shell using the LDA+U 
method. 

In the present paper we investigate the electronic struc- 
ture of the low-temperature phase using the LDA+J7 
approach in the tight-binding linear muffin-tin orbital 
(TB-LMTO) calculation scheme.SL22iSa The LDA+JJ 
calculations result in a charge and orbitally ordered 
insulator with an energy gap of 0.29 eV, which is in 
a good agreement with an experimental gap value of 
0.25 eV. From our results, we propose an orbital order 
parameter, defined as the difference between t2g major- 
ity/minority spin occupancies of Ti(l)'^+/Ti(3)'^"'" and 
Ti(2)'*+/Ti(4)^+ cations, respectively. This order param- 
eter is found to be quite large, although the total 3d 
charge difference between 3-1- and 4-1- cations, remains 
small. Also it is interesting to note that the total charge 
separation between 3-f and 4-1- Ti cations completely is 
lost due to efficient screening by the rearrangement of the 
other Ti electrons. In addition, we find a strong antifer- 
romagnetic coupling of J « 1700 K of the local moments 
within the dimcrizcd Ti'^+-Ti'^+ pairs, whereas an inter- 
pair coupling is only of « 40 K. This is in a good agree- 
ment with small and temperature independent magnetic 
susceptibility in the low-temperature phase of Ti407. 



II. COMPUTATIONAL DETAILS 

The present band-structure calculations have been per- 
formed for the low-temperature triclinic structure of 
Ti407.-^ The PI unit cell used in the calculations was 
constructed from the translation vectors of the original 
II ceU with a = 5.626 A, b = 7.202 A, c = 20.2608 
A, a = 67.90°, (3 = 57.69°, and 7 = 109.68° found 
at 115 K. The radii of muffin-tin spheres were taken 
as i?Tii-4 = 2.27 a.u., i?oi,03,04-6 = l-^S a.u., and 
-Ro2,07 = 1-66 a.u. Fifteen kinds of empty spheres were 
introduced to fill up the inter-atomic space. For simplic- 
ity we neglect small spin-orbit coupling and consider only 
a collinear spin case. 
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III. LSDA BAND STRUCTURE 

Our LSDA band structure calculations for the low- 
temperature PI structure confirmed the results of 
the previous worki2& The LSDA gives a nonmagnetic 
metallic solution with substantial charge separation be- 
tween crystallographically independent Ti(l)/Ti(3) and 
Ti(2)/Ti(4) cations. The lower part of the valence band 
below -3 eV is predominantly formed by O 2p states with 
a bonding hybridization with Ti 3d states. Crystal field 
splitting of the latter is roughly of 2.5 eV. Ti t2g states 
form the group of bands at and up to 2 eV above the 
Fermi energy whereas Ti eg states give a predominant 
contribution to the bands between 2.5 and 4.5 eV. Within 
the t2g group of bands the symmetry inequivalence of 
Ti(l)/Ti(3) and Ti(2)/Ti(4) sites leads to substantial 
t2g charge separation between these two groups of Ti 
atoms. In addition, an analysis of the partial density 
of states reveals significant bonding-antibonding splitting 
of dxy (in local cubic frame) states of about 1.5 eV for 
Ti(l)/Ti(3) cations, whereas Ti(2)/Ti(4) cations show a 
relatively weak substructure. This substantial bonding- 
antibonding splitting of Ti(l)/Ti(3) t2g states agrees 
well with the concept of formation of Ti'^+-Ti^+ spin- 
singlet pairs proposed earlier by MarezioiSLSS However, 
the LSDA calculations fail to reproduce an insulating 
spin-singlet ground state of the low-temperature phase 
of Ti407. Apparently, the electron-electron correlations, 
mainly in the Zd shell of Ti cations, play a significant 
role. 



IV. LDA+C/ RESULTS AND CHARGE 
ORDERING 

In order to take into account strong electronic cor- 
relations in Ti id shell we perform LDA-I-J7 calcula- 
tions for Ti407 in the low-temperature PI structure. In 
our calculations we use Coulomb intereaction parameter 
U — 3.0 eV and exchange coupling J = 0.8 eV taken in 
agreement with previous constrained LDA calculations.'^^ 
The LDA-|-[/ calculations result in a charge and orbitally 
ordered insulator with an energy gap of 0.29 eV (see 
Fig. 121) . This is in a strong contrast to the metallic so- 
lution with a substantial charge disproportionation be- 
tween crystallographically inequivalent Ti(l)/Ti(3) and 
Ti(2)/Ti(4) cations obtained by LSDA and in a reason- 
ably good agreement with an experimental gap value of 
0.25 eV. Note, however, that the charge and orbital or- 
der pattern remains exactly the same for U in the range 
2.5-4.5 eV, whereas the energy gap increases considerably 
up to 1.12 eV for C/=4.5eV. This remarkable increase of 
the gap value is accompanied by the enhancement of the 
spin magnetic moment from 0.56 up to 0.8 per 3-1- 
Ti(l)/Ti(3) cation as U is increased from 2.5 to 4.5 eV. 

In addition, we perform LDA-|-t/ calculations for high- 
temperature metallic phase of Ti407. In particular for U 
of 2.5eV a metallic self-consistent solution with substan- 
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FIG. 2: The total DOS obtained from LDA+(7 calculations 
with f/=3.0eV and J=0.8eV for the low-temperature Pi 
phase of Ti407. The top of the valence band is shown by 
dotted hnes. 

tial density of states (76 states/Ry) at the Fermi level has 
been found, whereas for J7 of 3 eV the LDA+J7 solution 
becomes unstable but remains metallic. With further in- 
crease of the U value the metallic solution collapses into 
insulating one. 

After self-consistency was achieved four crystallo- 
graphically independent Ti atoms are split out in two 
subgroups in respect to the spin magnetic moment per 
Ti site: Ti(l)/Ti(3) with a moment of 0.66/-0.67 ^is, re- 
spectively, and Ti(2)/Ti(4) with 0.04/-0.02 /i^. Thus, 
one of t2g majority /minority spin states of Ti(l)/Ti(3) 
becomes occupied (d^), whereas all other t2g states are 
pushed by strong Coulomb interaction above the Fermi 
level. In contrast, all t2g states of Ti(2) and Ti(4) are 
almost depopulated [d^) and form bands up to 2.5 eV 
above the Fermi level. The occupied Ti(l)/Ti(3) states 
are strongly localized and form a prominent structure 
with a band width of 0.25 eV just below the Fermi level 
(see Fig. ^ . The strong Coulomb interaction does not 
affect much the empty Ti Cg states, which give predomi- 
nant contribution between 2.5 and 4.5 eV. The obtained 
magnetic structure is almost antiferromagnetic with the 
spin magnetic moments within Ti(l)'^+-Ti(3)'^^ as well 
as Ti(2)''+-Ti(4)''+ pairs being of the same magnitude 
with opposite sign. 

An analysis of occupation matrices of Ti(l)'^+/Ti(3)''+ 
majority/minority 3rf spin states confirms substantial 
charge disproportionation within the Ti 3d shell. As 
shown in TableHl one of the t2g states of Ti^^^ cations [d^) 
is occupied with the occupation number of 0.74, whereas 
the remaining two t2g orbitals have a significantly smaller 
population of about 0.08. Thus, according to Ref. |^ we 
define an orbital order parameter as the largest difference 
between 3-1- and 4-1- Ti t2g populations which amounts 
to 66% of ideal ionic charge ordering model. The orbital 
order parameter clearly shows the existence of substan- 
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FIG. 3: The partial DOS for Ti(l)^+ and Ti(2)''+ cations are 
shown. The gap value of 0.29 eV was obtained by LDA+(/ 
with U—3.0eV and J=0.8eV. The Fermi level is shown by 
dotted line. 



tial charge disproportionation in the Ti 3c? shell of Ti407 
which is remarkable because of the complete lack of the 
total charge separation (see column q in TableQJ between 
3+ and 4+ Ti cations. The occupation matrices analy- 
sis shows that the change of the t2g occupations is very 
efficiently screened by the rearrangement of the other Ti 
electrons. A significant portion of the screening charge 
is provided by Ti eg states due to formation of relatively 
strong a bonds with O 2p states, which results in appre- 
ciable contribution of the former to the occupied part of 
the valence band. Ti 4s and Ap states give additional 
contributions to the screening of the difference in t2g oc- 
cupations which leads to complete loss of the dispropor- 
tionation between the charges at 3-1- and 4-1- Ti sites. 

The occupied t2g Ti'^^ states are predominantly of dxy 
character in the local cubic frame (according to that we 
later mark the orbital as dxy orbital). This is illustrated 
in Fig.0I which shows the angular distribution of the ma- 
jority and minority spin 3c? electron density of Ti cations, 
marked by red and cyan color (or light and dark on the 
black and white image), respectively!?^ Since Ti(l)'^+ and 
Ti(3)'^"*' cations are antiferromagnetically coupled, the 
obtained ferro-orbital order is consistent with the forma- 
tion of a bonding spin-singlet state from the d^y orbitals 
of two neighboring Ti(l) and Ti(3) sites. The orientation 



TABLE L Total (g) and i-projected {qs.p,d) charges, mag- 
netic moments (M), and occupation of the most populated 
t2g orbitals (n) calculated for inequivalent Ti atoms in the 
low-temperature PI phase of Ti407. 



Ti ion 


<? 


qs 




Qd 


M (^ib) 


t2g orbital 


n 


Ti(l)^+ 
Ti(2)''+ 


2.27 


0.18 


0.27 


1.83 


0.66 


dxyl 


0.74 


2.22 


0.22 


0.33 


1.68 


0.04 




0.08 


Ti(3)='+ 
Ti(4)''+ 


2.16 


0.18 


0.25 


1.74 


-0.67 


dxy I 


0.73 


2.16 


0.21 


0.33 


1.62 


-0.02 




0.07 




Ti 1 



Ti 3 



FIG. 4: (color online) Structure of Ti407 showing the angu- 
lar distribution of the majority and minority spin 3c! electron 
density of Ti cations. Red and cyan (light and dark, respec- 
tively, on the black and white image) orbitals correspond to 
the majority and minority 3c? spin states, respectively. Oxy- 
gen atoms are shown by small spheres. The size of orbital 
corresponds to its occupancy. X-Y-Z coordinate system cor- 
responds to the local cubic frame. 



of occupied Ti'^+ t2g orbitals is consistent with the largest 
average Ti-0 distance in the plane of t2g orbitals. As 
shown in TableHIlthe average Ti(l)-0 distance (2.061 A) 
in the plane of d^y orbital is considerably larger than av- 
erage distances in the other two yz and zx planes (2.032 
and 2.045 A, respectively). The same is also true for the 
Ti(3) cation but in this case the variation of the average 
Ti(3)-0 distances is much smaller (2.047 vs 2.041 and 
2.042 A) and, as a consequence, the out-of-plane rota- 
tion of the occupied t2g minority spin orbital is stronger. 

In addition, hopping matrix elements were evaluated 
via Fourier transformation from reciprocal to real space 
of the Ti t2g LDA Wannier Hamiltonian.'^^ Remarkably, 
for the low-temperature phase the Ti(l)-Ti(3) intra-pair 
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TABLE II: The averaged Ti-0 distances in the plane of t2g or- 
bitals (dorb.) and in the oxygen octahedra (dav.) for PI struc- 
ture of Ti407. The occupied orbital of the Ti(l) and 
Ti(3) 3+ cations is predominantly of d^y character. 



Ti atom 


orbital 


darh. (A) 


da„. (A) 


Ti(l) 


dxy 


2.061 


2.046 




dyz 


2.032 






dzx 


2.045 




Ti(2) 


dxy 


2.012 


2.000 




dyz 


1.976 






dzx 


2.013 




Ti(3) 


dxy 


2.047 


2.043 




dyz 


2.041 






dzx 


2.042 




Ti(4) 


dxy 


1.973 


1.977 




dyz 


1.976 






dzx 


1.981 





dxy-dxy hopping matrix element is found to be of 0.61 
eV, whereas all other hoppings are 3-4 times smaller. 
This strong inhomogeneity of the hopping matrix ele- 
ments disappears in the high- temperature phase. Thus, 
according to our calculations hopping elements in the 
high-temperature phase are 0.23, 0.21, 0.39, and 0.33 eV 
between 1-3, 2-4, 3-3, and 4-4 Ti sites, respectively. 

Estimation of exchange interaction parameters via the 
variation of the ground state energy with respect to the 
magnetic moment rotation angl o^^i^^ results in a strong 
antiferromagnetic coupling of -1696 K between Ti(l)'^+ 
and Ti(3)'^"'" cationsi^ All other couplings are two orders 
of magnitudes smaller. This indicates a possible forma- 
tion of the spin-singlet pairs via direct antiferromagnetic 
exchange between neighboring Ti(l) and Ti(3) sites. The 
contribution of the superexchange via O p orbitals to the 
Ti(l)-Ti(3) exchange coupling is found to be negligible. 
This was verified by calculating the exchange coupling 
constants with the sub-blocks of the LMTO Hamiltonian 
responsible for the Ti-0 hybridization being set to zero. 



This calculation gave qualitatively same results for the 
exchange constants although the possibility for the su- 
perexchange via O p orbitals was eliminated. 



V. SUMMARY AND CONCLUSIONS 

In summary, in the present LDA+J7 study of the low- 
temperature P\ phase of Ti407 wc found a charge or- 
dered insulating solution with an energy gap of 0.29 eV. 
The total ?>d charge separation is small (less than 0.14), 
whereas the orbital order parameter defined as the differ- 
ence between t2g occupancies of Ti^+ and Ti^'^ cations 
is large and gives direct evidence for charge ordering. Ti 
4s and Ap states give a strong contribution to the static 
"screening" of the total 3£i charge separation. This ef- 
fective charge screening leads to complete loss of the dis- 
proportionation between the charges at 3-1- and A+ Ti 
sites. The occupied t2g states of Ti^+ cations are pre- 
dominantly of dxy character (in the local cubic frame) 
and form a spin-singlet molecular orbital via strong di- 
rect antiferromagnetic exchange coupling between neigh- 
boring Ti(l) and Ti(3) sites of J w 1700 K, whereas the 
role of superexchange is found to be negligible. This is 
in a good agreement with small and temperature inde- 
pendent magnetic susceptibility in the low-temperature 
phase of Ti407. 
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The exchange couphng parameter J represents the effective 
pair excfiange interaction between Ti atoms with effective 
Heisenberg Hamiltonian H = — X)i>j Jij^i ■ Sj, where Si 
and Sj are spin magnetic moment vectors at site i and j. 
Positive (negative) values of J correspond to the ferromag- 
netic (antiferromagnetic) coupling between sites. 



